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[57] ABSTRACT 

A frequency-domain block-adaptive digital filter 
(FDAF) having a finite impulse response of length N 
for filtering a time-domain input signal in accordance 
with the overlap-save method includes window means 
(11) for obtaining modifications (B(p;m)) of the 2N 
frequency-domain weighting factors (W(p;m)) from 
correlation products (A(p;m)). A known FDAF of this 
type contains five 2N-point FFT's, two of which are 
used in the window means (11). By utilizing a special 
time-domain window function which can be imple- 
mented very efficiently in the window means (11) with 
the aid of a frequency-domain convolution, a FDAF of 
this type containing only three 2N-point FFT's is ob- 
tained whose convergence properties are comparable to 
those of the known FDAF" containing five 2N-point 
FFT's. 

1 Claim, 2 Drawing Sheets 




05/14/2004, EAST Version: 1.4.1 



U.S. Patent Feb. 21, 1989 Sheet 1 of 2 4,807,173 



✓Sectioning 

x(k))^1 x li.m)^ X(p ; m) 



on 

2N 



CONJUGATION 



H H 



^o/|x<p;m)|< 
2p(p ; m ' 



U 13 W(p ; m) 



Y(P;(Tj| 




Alp ; m 



I0FT 


\ . 


2N 





T 
17 



n 

o(i;m) . bli.m) Blp.m) 



• 


\ . 


OFI 






2N 



IB 

g(i) 



A 



A(p ; m) B(p : m ) 



19 



3 



20 



FIG.2Q 



G(p) 

FIG.2b 



05/14/2004, EAST Version: 1.4.1 



US. Patent Feb. 21,1989 Sheet 2 of 2 4,807,173 




05/14/2004, EAST Version: 1.4.1 



4,807,173 

1 2 

_ „ nals in the digital filter are represented. With the cus- 

FREQUENCY.D^^ tomary( valid ^Jptvns about the statistic 

DIGITAL FILTER independence of the different quantities in the filter, the 

BACKGROUND OF thp rwvp wrrnw c we « htta * fac ! ors ^ have to same variances when no 
BACKGROUND OF THE INVENTION 5 use is made of window functions in the adaptation loop 

The invention relates to a frequency-domain block- of to This implies that at the same convergence 
adaptive digital filter having a finite impulse response of rate of the adaptive filter (that is to say the same gain 

length N for filtering a time-domain input signal in factor in the adaption algorithm) using 2N instead of N 

accordance with the overlap-save method. weighting factors results in an increase of the final mis- 

A frequency-domain adaptive filter (FDAF) having 10 alignment noise factor of the filter by 3 dB f since the 

such a structure is disclosed in the article "A Unified final misalignment noise factor is determined by the sum 

Approach to Time- and Frequency-Domain Realization of the variances of the weighting factors In practice, 

™^ daptive Digital Rtcrs " by G ' A " aark et d - m the «™ factor m adaptation algorithm is chosen 
™ J™*- Acoust, Speech, Signal Processing, Vol. such that a predetermined value of the final misalign- 
ASSP-31, No. 5, October 1983, pages 1073-1083, more «■ ment noise factor is not exceeded. In order to compen- 
I w r u a * ... for to increase of the final misalignment noise 

In the field of speech and data transmission, tune- factor in an FDAF, this gain factor must be halved, 

£"TJ ? aP * ^ " «d which ~ to convergence rate also to be halved 

TOAFW^ltl^H ^ hcati0 ^ hese m whereas in the majority of applications a highest possi: 

TOAP s are implemented as adaptive transversal filters, 20 w ^ fate . 

m which a "least-mean-square" (IMS) algorithm is Q.M.Ak,,ri.A H T u . , 
usedforadapthgtheweigha^ pro^m^^ 

impulse response assumes large values, as is the case / • * u 3 : the r m ° dlflca ' 

withapp^^ S/de^ 

mented as a transversal filter has the disadvantage that 25 bem S directly from the second multiplier 

the complexity in terms of arithmetic operations (multi- meanS /- but b J "P 8 T d ° W meanS for P erformin g ™ 
plying and adding) per output sample increases linearly °P er * Uon W h ? e " domam equivalent is a multiplica- 
with the filter length N. In addition, the TD AF imple- by a rec ^ n Sy lar win <*ow function of length 2N 

mented as a transversal filter has a low convergence w , constrams the last N components to be zero. An 
rate for highly correlated input signals such as speech 30 ™Plementation of this window function in the time- 
and certain types of data, since the convergence rate domain requires the use of 2 DFT's, namely an inverse 
decreases with an increasing ratio of the maximum to *° r to transformation to the time-domain and a 

minimum eigenvalues of the correlation matrix of the DFT to effect tne transformation to the frequency- 
input signal (see, for example, C. W. K. Gritton and D. domain after multiplication by the time-domain window 
W. Lin, "Echo Cancellation Algorithms", IEEE ASSP 35 function. An alternative implementation is based on the 
Magazine, April 1984, pp. 30-38, in particular pp. consideration that a multiplication in the time-domain is 
. 32/33). equivalent to a convolution in the frequency-domain 

The use of frequency-domain adaptive filters to components of the DFT of length 2N of this 

(FDAF) provides the possibility to significantly im- time-domain window function. For high values N this 
prove ^he convergence properties for highly correlated 40 alternative implementation in the frequency-domain is 
input signals, as for any of the substantially orthogonal not attractive, as its computational complexity per com- 
frequency-domain components of the input signal the ponent increases linearly with N, whereas in the first- 
gain factor in the adaptation-algorithm can be normal- mentioned time-domain implementation this complexity 
ized in a simple way in accordance with the power of becomes proportional to the logarithm of N when the 2 
the relevant frequency component. For the most effi- 45 DFT's are implemented as FFTs. Thus, a preferred 
cient implementation of a FDAF having an impulse implementation of the known solution results in an 
response of length N, use is made of Discrete Fourier FDAF containing a total of 5 DFTs implemented as 
Transforms (DFT) of length 2N of 2N weighting fac- FFT's. 
tors to ensure that circular convolutions and correla- 
tions, computed with the aid of DFT's, are equivalent 50 SUMMARY OF THE INVENTION 
to the desired linear convolutions and correlations The invention has for its object to provide a frequen- 
when the sectioning method is performed correctly. cy-domain block-adaptive digital filter of the type de- 
For large values N the computational complexity can, fined in the opening paragraph of section (A) whose 
however, be significantly reduced in terms of arithmetic computational complexity is significantly reduced com- 
operations per output sample by utilizing efficient im- 55 pared to that of the known FDAF containing 5 DFT's, 
Tte^S^* 118 DFT u kn ° w ?. » T 351 Fourier by utilizing a priori information about the global shape 
Transform (FFT), as a result of which this complexity 0 f the impulse response to be modelled, but whose con- 
^ Proportional to the logarithm of the filter vergence behaviour is comparable to that of the known 

tt? i. ™ AC _ , . FDAF containing 5 DFT's. 

There where a TD AF needs only have N weighting 60 The frequency-domain block-adapti ve digital filter 
■|« ^ap™ m P^. res P°" se N, the equiva- according to the invention is characterized, fa that the 

other types of signals superimposed on the reference ^ of len ^ h 2 N defined hv nmcuon 
signal and because of the precision (that is to say word gW ^ V * 
length or number of bits) with which the different sig- . g(k)Mi) [i+co* {(k~k 0 yn/N)) 
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for k=0, 1, • • . , 2N— 1, whereko is a constant with 
0^ko<N 

SHORT DESCRIPTION OF THE DRAWINGS , 

The invention and its advantages will now be de- 
scribed in greater detail by way of example with refer- 
ence to the accompanying drawings. Therein: 

FIG. 1 shows the general block diagram of an FDAF 
including window means and utilizing the overlap-save 10 
method; 

FIGS. 2a and 2b show the block diagram of the win- 
dow means for the implementation of a window func- 
tion in the time-domain and the implementation, equiva- 
lent thereto, in the frequency-domain, respectively; 15 

FIGS. 3a and 36 show the time-domain window func- 
tion according to the afore-mentioned prior art and the 
invention, respectively; 

FIG. 4 is a time diagram to illustrate the convergence 
behaviour of an FDAF in accordance with the afore- 20 
mentioned prior art and the invention, respectively. 

DESCRIPTION OF THE EMBODIMENTS 

FIG. 1 shows the general block diagram of an FDAF 
having a finite impulse response of length N for filtering 25 
a time-domain digital input signal x(k) in accordance 
with the overlap-save method. Double-line signal paths 
in FIG. 1 indicate paths in the frequency-domain and 
single-line signal paths indicate paths in the time- 



signal block m. With the aid of transformation means 4 
for. performing an inverse 2N-point DFT each output 
signal block is transformed to the time-domain;, the 
resultant time-domain points of block m being denoted 
shown as y(i,*m) with i=0, 1, . . . 2N— 1. Since the 
weighting factors W(p;m) may be considered as being 
points of a 2N-point DFT performed on time-domain 
weighting factors w(i,*m) which represent values of the 
impulse response w(k) during block m, the multiplica- 
tion in multiplier means 3 in the time-domain corre- 
sponds to the circular convolution of the input signal 
i(k) during block m with the impulse response w(k) 
during block m. The desired time-domain output signal 
y(k) is, however, the linear convolution of x(k) with 
w(k). In accordance with the overlap-save method, this 
desired output signal y(k) is obtained by applying the 
time-domain components y(i;m) of this circular convo- 
lution for each block to sectioning means 5 in which, 
using parallel-to-serial conversion, the first N point 
y(i;m) with i=0, 1, . . . , N— 1 are discarded and the last 
N points y(i;m) with i=N, N+ 1, . . . , 2N- 1 are passed 
on as the output signal y(k), as is symbolically shown in 
FIG. 1. 

For the adaptation of the frequency-domain 
weighting factors W(p;m) on a block-by-block basis, use 
is made of a "least mean-square" (LMS) algorithm. In 
accordance with this algorithm, these weighting factors 
W(p;m) are modified as long as there is correlation 
between the input signal x(k) and an error signal r(k) 



domain. Transformation from the time-domain to the 30 given by the difference between the reference signal 



frequency-domain and vice-versa is effected with the 
aid of the Discrete Fourier transform (DFT) and the 
Inverse Discrete Fourier transform (IDFT), respec- 
tively, both having a length 2N. In literature these 
transforms are known as 2N-point DFTs, wherein 35 
"point" may refer to both a discrete time-domain com- 
ponent and to a discrete frequency-domain component. 
To differentiate between time-domain and frequency- 
domain signals, time-domain signals are written in low- 
er-case letters and frequency-domain signals in upper- 40 
case letters. 

The FDAF shown in FIG. 1 has for its object to 
derive, at discrete instants k, from the input signal x(k) 
an output signal y(k) which equals a reference signal 
d(k) as well as possible. In many case, for example, 45 
when the FDAF forms part of an echo cancellor, this 
reference signal d(k) may be assumed to be the linear 
convolution of the input signal x(k) with ah impulse 
response h(k) of length N whose shape is not accurately 



d(k) and the output signal y(k). This differential signal 
r(k)=d(k)— y(k) is obtained with the aid of an adder 6. 
The overlap-save method for determining this correla- 
tion between the signals x(k) and r(k) implies that in 
FIG. 1 the error signal r(k) is applied to sectioning 
means 7 in order to be segmented into blocks of length 
2N with the aid of serial-to-parallel conversion, each 
block overlapping its preceding block over a length N 
and the forst portion of length N of each block being 
constrained to be zero, as is symbolically shown in FIG. 
1. The points of error signal block m are denoted as 
r(i;m) with i=0, 1, . . . , 2N— 1, it holding for i=0 t 1, ♦ 
. . , N— 1 that r(i;m)=0. Using transformation means 8 
for performing a 2N-point DFT each error signal block 
is transformed to the frequency-domain, the frequency- 
domain points of block m thus obtained being denoted 
as R(p;m) with p=0, 1, . . . , 2N— I. In addition, the 
frequency-domain components X(p;m) of input signal 
block m are applied to conjugation means 9 for produc- 



known. Then the FDAF has for its task to make its 50 ing the complex conjugate value X* (p;m) of each corn- 



impulse response w(k) equal to this impulse response 
h(k) as well as possible. 

To that end, in FIG. 1 the input signal x(k) is applied 
to sectioning means 1 in order to be segmented into 
blocks of length 2N with the aid of serial-to-parallel 55 
conversion, each block overlapping its preceding block 
over a length N, as is shown symbolically in FIG. 1 
The points of the input signal block having block num- 
ber m are denoted as x(i;m) with i=0, 1, . . . , 2N— 1. 
With the aid of transformation means 2 for performing 60 
a 2N-point DFT each input signal block is transformed 
to the frequency-domain, the frequency-domain points 
of block m thus obtained being denoted as X(p;m) with 
p=0, 1 . . . , 2N— 1. In multiplier means 3 each frequen- 
cy-dpmain component X(p;m) is multiplied with an 65 
associated frequency-domain weighting factor W(p;m) 
for forming products X(p;m) W(p;m) representing the 
frequency-domain components Y(p;m) of the .output 



ponent X(p;m). Each conjugated component X* (p;m) 
is multiplied in multiplier means 10 with the associated 
component R(p;m) for forming frequency domain prod- 
ucts X* (p;m) R(p;m) which correspond to the time- 
domain correlation between the input signal x(k) and 
the error signal r(k) during block m. In addition, each 
product X* (p;m) R(p;m) is multiplied by an amount 
2fi(p;m), where jx(p;m) is the gain factor in the adapta- 
tion algorithm, as a result of which a product A(p;m);32 
2fi(p;m) X*(p;m) R(p;m) is formed which determines 
the modification of the frequency-domain weighting 
factor W(p;m). In the FDAF shown in FIG. 1, these 
products A(p;m) are applied to window means 11 for 
obtaining the ultimate modifications 8(p;m) of the 
weighting factors W(p;m) which are formed by means 
of accumulator means 12. These accumulator means 12 
include a memory 13 for storing the weighting factors 
W(p;m) of block m and an adder 14 for forming the sum 
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For performing its main function (effecting that its 5 I ♦ ^ ^ f onv ^ nc f 15 ^wn as final misalign- 
impulse response w^n™ to ^SuS £ 

sponse h(k) to be modelled, a^FD^ uSf t£ determm ^ by the sum of the weighting factor van- 
overbp-sav'e method need not to u£ SfSlS 5~£ t^SH? ? ,? ^T^*^ 
means 11 of FIG. 1. In that case the modifications * j? fa f °V ? T ,°[ th JL 2N fc^ 0 ™* 
B(p;m) of the weighting factors W(p;m) are given by io w k m > °l *** roAF are actua "y 

the products A(p;m) and the adaptatSn algorithm may m P he * that the final value fi for the FD AF 

be written as: 13 mdeed unnecessarily higher by a factor of two (3 dB) 

than the final value 0 for the equivalent TDAF having 
wfom+ 1)= i^™)+2p<p;m)x*<p;m)R(p;m\ (i) om *y N time-domain weighting factors. For the custom- 

ary values of the adaptation factor a this final value £ is 
When the input signals are not or only weakly corre- substantially proportional to a and the same holds for 
lated, the gain factor ji(p;m) for each weighting factor convergence rate of 0(m) as long as £(m) is much 

W(p;m) may have the same constant value a which is greater than £. In practice, the adaptation factor a is 
independent of the block number m (this constant a is chosen such that a predetermined final value 0 is not 
known as the adaptation factor of the algorithm). For - 0 exceeded. This implies that the adaptation factor a for 
highly correlated mput signals the convergence rate can the FDAF must actually be halved unnecessarily to 
be significantly increased by decorrelating these input satisfy this prescription resulting in the convergence 
signals, which can be effected by normalizing their rate also being halved, which is undesirable 
power spectrum (see for example, the aforementioned This problem can be solved by including the window 
article by Gntton and Lm page 36). Since the frequen- means 11 in the FDAF of FIG. 1. In accordance with 
S^SES 1 S** 000 ? X <P ;m > *? *L™*y mailable in 25 the article by Clark et al. mentioned in the foregoing, 
an FDAF, the normalization can be effected in a simple these window means 11 can be arranged for performing 

£S^«ff ° f \ hC P ° Wer |X(p;m)l ° f s*™^ b ? ^tiplier means 10 that h eqdvaleitTo 

the relevant component X(p;m): ^ the J ociated domain prXte a^rS 

tfm)=*/\Xtom)p m 3 r t Ctangular willdow Unction of length 2N which 

w forces the last N time-domain products aftm) to be 

This possibility is shown in FIG. 1 by including means FIGS * 20 *?* 2b S , how the block dia S ram of these 

IS in which for each applied component Xfam) the ^ 11 ™° h M th * y m Ascribed in the 
righthand side of formula (2) is detennined which, op. M T^r^K* (see nG * 3) ' In the window means 

tionally after smoothing on a block-by-block basis, with 2 \? -T Ume - domain window function g(k) is 

the aid of a simple recursive filter, is used as a gain rea ^ zed v* 1 ™ tne of transformation means 17 for 

factor /t(p;m) for multiplying by 2/i(p;m) in the multi- JS™ 01118 m inverse 2N *P oint DFT converting the 

plier means 10. 2N frequency-domain products A(p;m) into the 2N 
As has already been mentioned in the foregoing, an 40 ^o^ated time-domain products a(i;m), a multiplier 18 

FDAF must utilize 2N-point DFTs and 2N frequency- mu "iplying each time-domain product a(i;ro) with the 

domain weighting factors W(p;m) for providing an Y alue s ® of window function g(k) for k=i to obtain 2N 

impulse response w(k) of a length N. These weighting time-domain products b(imi)=g(i) a(i;ra) and transfor- 

factors W(p;m) may be considered to be the points of a mat ion means 19 for performing a 2N-point DFT con- 
2N-point DFT performed on 2N time-domain 45 verting the 2N time-domain products b(i;m) into the 

weighting factors w&m). This creates a problem be- associated 2N frequency-domain products B(p;m). 

cause the FDAF has only N time-domain weighting Tllese P roduct s B(p;m) constitute the modifications of 

factors w(i;m) which represent the N values w(k) with the frequency-domain weighting factors W(p;m) 

k=0, 1, . . . , N-l of the impulse response of the applied to the accumulator means 12 of FIG. 1. The 
FDAF, so that the remaining N time-domain weighting 50 window means 11 of FIG. 2b are based on the consider- 

factors W&m) with i;32 N, N+l. . . . , 2N-1 are actu- ation H»t multiplying by g(i) in the time-domain as 

ally superfluous. It is likewise described that, after con- performed in FIG. 2a is equivalent to a convolution in 

vergence, the 2N weighting factors W(p;m) and conse- the frequency-domain and consequently consist of con- 

quently also the 2N weighting factors w(i;m) will con- volution means 20 for performing the (circular) convo- 

tuiue to fluctuate around their final values as a result of 55 wtion of the frequency-domain products A(p;m) with 

ttus finite precision of the signal representation in the the frequency-domain components G(p) obtained by 

FDAF and because of the presence of noise and any effecting a 2N-point DFT on the 2N points g(i) of the 

other types of signals superimposed on the reference time-domain window function g(k) for k^i 

agr^d(k) This superpositioning is symbolically shown The rectangular window function g(k) of length 2N 

in FIG. 1 by a broken-line adder 16 inserted in the path 60 as described in the article by Clark et al. is eiven bv the 

of the reference signal d(k) and receiving a signal s(k) formula: 
representing noise and any other type of signals. The 
variances of the fluctuating weighting factors will all 

have the same value when the FDAF does not include rfJt * _ ( l * = o, l n - i (3) 

window means 11 and the gain factors have been 65 \o k = N,N+i w- i 

chosen in the manner described in the preceding para- 

graph. An important parameter for the convergence and is shown in FIG. 3a For the sake of simolicitv the 

behaviour of a block-adaptive filter (FDAF and urftacnceofto 
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gence behaviour of the FDAF will now be described in higher value than for the case in which the FDAF does 
greater detail for the case in which the window means not include window means 11. When this twice higher 
11 of FIG. la are used in FIG. 1 and the memory 13 has value a =2- 5 is also used in the FDAF without win- 
zero content at the beginning of the adaptation. This last dow means 11, )3(m) varies in accordance with curve c 
fact implies that for m=0 the 2N frequency-domain 5 in FIG. 4 which coincides with curve b until (m) has 
weighting factors W(p;m) satisfy W(p;0)=0 and conse- reached the value £(m)= — 15 dB, which at this value 
quently the associated time-domain weighting factors of a constitutes the final value ft consequently this final 
w(i^n) also satisfy w(i0)«=0. each of the 2N weighting value /3= — 15 dB exceeds the prescribed final value 
factors w(i;m) for m^O may be considered to be the /3= — 18 dB by 3 dB. 

result of the block-by-block accumulation of its time- 10 As regards the practical implementation of window 
domain modifications starting from m=0 and when the means 11 for providing window function g(k) as defined 
window means 11 of FIG. 2a are used, these time- by the formula (3), the embodiment of FIG. 2a should 
domain modifications are formed by the products be preferred over the embodiment of FIG. 2b. more 
b(i;m)— g(i)a(i;m) at the output of multiplier 18. By specifically for large values N, as the computational 
using window function g(k) as defined by formula (3), 15 complexity, expressed in the number of real multiplica- 
the modifications b(i;m) with i=N, N+ I, . . . » 2N— 1 tions and additions required for obtaining the 2N modi- 
and m=0, .1, 2, ... are constrained to be zero and conse- fications B(p;m) of the weighting factors W(p;m), is in 
quently the weighting factors w(i;m) with i-N, N+l, the case of FIG, 2a of a lower order than for the case 
... , , 2N— 1 for m^O are also constrained to be zero illustrated in FIG. 2b, FIG. 2a requires two 2N-point 
because then w(i;m)=w(i;0) for these values i and be- 20 DFTs which can efficiently be implemented as 2N- 
cause it holds for all values i that w(i;0)=0. This implies point FFTs, so that for large values N .the number of 
that of the 2N weighting factors wftm) of the FDAF arithmetic operations is of the order N log(N), the mul- 
the N actually superfluous weighting factors w(i;m) tiplication by g(k) in multiplier 18 not contributing to 
with i=N, N+ 1, . . . , 2N- 1 will not fluctuate around this number of operations since $j(k) only assumes the 
the constrained value zero and consequently will not 25 values I and 0. In FIG. 2b the 2N values A(p;m) are 
contribute to the final value £ of the parameter /3(m). convolved with the values G(p) obtained by performing 
Halving the adaptation factor a, necessary in the a 2N-point DFT on the 2N points g(i) of g(k) in accor- 
FDAF without window means 11 to prevent a pre- dance with formula (3) for k=i. As can be easily 
scribed final value & from being exceeded, needs conse- checked, G(p) has the value G(p)=0 for p=2, 4, 6, . . . 
quently not to be effected in the FDAF including win- 30 ^ a value G(p)#0 for p =0 and p = 1, 3, 5, ... So FIG. 
dow means 11 for realizing window function g(k) in 2b requires the convolution of 2N values A(p;m) with 
accordance with formula (3) and consequently there is (NOl) values G(p), so that for large values N the num- 
also no need to halve the convergence rate. ber of arithmetic operations is of the order of N 2 and 

The convergence behaviour of the described consequently of a higher order than for the case illus- 
FDAFs will now be illustrated with reference to simu- 35 trated in FIG. 2a. 

lation results for the case in which the impulse response The FDAF's described so far do not utilize a priori 
h(k) of length N=32 to be modelled has the global information about the usually well-known global shape 
shape of an exponentially decreasing function, the input 0 f the impulse response h(k) to be modelled and only 
signal x(k) is a pseudo-ternary data signal in accordance utilize the a priori information about the length N of this 
with the AMI-code and is, therefore, a highly corre- 40 impulse response. In contrast thereto, the FDAF ac- 
lated signal, the value )3=2-* (= -18 dB) is chosen as cording to the invention does indeed utilize a priori 
the prescribed final value of the parameter 0(m) and the information about the global shape of the impulse re- 
value £(0)=2 7 (= +21 dB) is chosen as the initial value. sponse h(k) to be modelled, as usually it is not only 
The simulation results are represented in a highly styl- known of this impulse response h(k) that the N values 
ized form in FIG. 4, which shows the parameter /5(m) as 45 h(k) for k=0, 1, , . , t N- 1 differ significantly from zero 
a function of the number m of the iterations of the adiap- and are substantially zero for k^N, but it is also 
tation algorithm. This stylizing relates to both the con- roughly known for which values k the largest ampli- 
vergence rate of /3(m) and also to the details of the tudes |h(k)l occur. In addition, it is known, that the 
variations of /J(m) for consecutive values of m. More amplitudes (h(k)| have roughly a decreasing character 
specifically, the convergence rate of fi(m) in FIG. 4 has 50 for values k which increasingly differ from the values k 
a constant value until the final value ft is reached and y/ith the largest amplitudes |h(k)|. The invention now 
thereafter has zero value, whereas this convergence utilizes this a priori information by arranging the win- 
rate is actually already noticeably smaller at values of dow means 11 for providing a time-domain window 
/5(m)=2- 2 (= -6 dB) than at the initial value fi(0) and function g(k) of a length 2N which is given by the for- 
still further decreases at still lower values of J3(m); and, 55 m ula: 
additionally, the details of the variation of 0(m) for 

consecutive values of m have been omitted in FIG. 4 to s(*)=(J)[l+cos ((A-fc^/N)] (4) 

prevent the noisy character of these details from ob- 
scuring the image of the global shape of £(m) to too for k-0, 1, . . . m 2N- 1, where k<, is a constant with 
high an extent. 60 0^ko<N. By choosing this constant ko in the range of 

Curve a in FIG. 4 relates to the FDAF without win- values k having the largest amplitudes |h(k)| it is ac- 
dow means 11, the adaptation factor having a value complished that the window function g(k) defined in 
a =2~ 6 for reaching the prescribed final value /? = — 18 formula (4) has the same global shape as the impulse 
dB. Curve b in FIG. 4 relates to the FDAF including response h(k) to be modelled. In addition, the window 
window means 11 for implementing window function 65 means 11 are arranged for implementing -the window 
g(k) as defined in formula (3), the adaptation factor function g(k) defined by formula (4) in the frequency- 
having a value a=2 _5 for reaching the prescribed final domain, that is to say for the convolution of the 2N 
value 0= — 18 dB; this adaptation factor a has a twice frequency-domain products A(p;m) produced by multi- 
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plier means 10 in FIG. 1 with the frequency-domain 
components G(p) obtained by effecting a 2N-point 
DFT on the 2N points g(i) of g(k) defined by formula 
(4) for k=i. As is easy to check, G(p) has the following 
values: 5 

G(0)=N 

G(»)=(W2)[cos(Aoir/A)-ysin<WA0] • 

<^-l)«(W2)[co9(^/iS0+yu(WA)] (5) 10 

G(p)=0, (5 ) 

Implementing time-domain window function g(k) de- 
fined in formula (4) in the frequency-domain with the 15 
aid of window means 11 as shown in FIG. 2b conse- 
quently requires the convolution of 2N values A(p;m) 
with only 3 values G(p), namely the real value G(0) and 
the mutually conjugate complex values G(l) and 
G(2N- 1), so that for high values N the number of 20 
arithmetic operations required for the 2N modifications 
B(p;m) of the weighting factors W(p;m) is only of the 
order N and consequently not of the order N 2 , as in the 
case of the known window function g(k) as defined in 
formula (3). Using the measures according to the inven- 25 
tion consequently accomplishes a significant reduction 
of the computational complexity. In this respect it 
should be noted that this reduction of the computational 
complexity is not achieved when the window function 
g(k) of formula (4) is implemented with the aid of win- 30 
dow means 11 as shown in FIG. 2a, because multiplying 
by g(k) then indeed contributes to the number of arith- 
metic operations, but particularly because the computa- 
tional co mplexi ty of the two 2N-point DFTs imple- 
mented as FFT does not change, so that for large values 35 
N the number of computational operations remains of 
the order N log(N). 

Now it will be explained that compared to the pre- 
ferred embodiment of the known FDAF with a total of 
5 2N-point DFTs implemented as FFT the significant 40 
reduction in the computational complexity of the* 
FDAF according to the invention is not accompanied 
be a deterioration in the convergence behaviour. For 
simplicity, this explanation will be given with reference 
to observations in the time-domain. 45 

In the foregoing it has already been stated that the 
FDAF actually has 2N time-domain weighting factors 

w(i;m) with i=0, 1 2N- 1. When the FDAF does 

not include window means 11 and the gain factors 
M(PJni) have been chosen in accordance with formula 50 
(2), the rate at which a weighting factor w(i;m) con- 
verges to the associated value h(i) of the impulse re- 
sponse h(k) to be modelled, will have for all weighting 
factors w(i;m) the same value which is proportional to 
the adaptation factor a. Because of the use of time- 55 
domain window functions g(k) the convergence rate of 
a weighting factor w&m) becomes proportional to g(i). 
The convergence rate of the FDAF is a function of the 
convergence rate of all the weighting factors w(i;m), 
but will predominantly be determined by the conver- 60 
gence rate of those weighting factors w(i;m) which 
have to converge to the highest values h(i) of the im- 
pulse response h(k) to be modelled. By now choosing 
the constant k<, in formula (4) in the range of the values 
k having the largest amplitudes |h(k)| it is accom- 65 
plished that the value g(k) in that range is substantially 
equal to 1, that is to say it has substantially the same 
value as in the case in which no window functions g(k) 
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arc used (in fact, this latter case can also be character- 
ized by a window function g(k)= 1 for k=0, 1, . . . , 
2N-1). Chossing the window function g(k) in accor- 
dance with formula (4) then results in a convergence 
rate of the FDAF of the invention substantially equal to 
that of the FDAF without window means 11 when the 
adaptation factor has the same value for both cases. 

As has also been mentioned in the foregoing, the 
variance of a weighting factor w(i;m) after convergence 
has the same value for all the weighting factors w(i;m) 
when the FDAF does not include window means and 
the gain factors fi(p;m) have been chosen in accordance 
with formula (2), For the usual values of the adaptation 
factor a the value of this variance of w(i;m) is substan- 
tially proportional to a. However, when the window 
function g(k) is used, the variance of a weighting factor 
w(i;m) after convergence becomes proportional to 
ag(i); The final value 0 of the parameter (m) is substan- 
tially determined by the sum of the variances of the 
weighting factors w(i;m) after convergence, so that 
when the window functions g(k) are used the final value 
0 is substantially proportional to the quantity 



2N-\ 

1 A M 



(6) 



for the customary values of the adaptation factor a. For 
the FDAF without window means 11 the quantity of 
formula (6) has the value 2aN and for the FDAF with 
window means 11 this quantity has the value N when 
the window function g(k) is in accordance with formula 
(3), but also when the window function g(k) is in accor- 
dance with formula (4). This can be checked in a simple 
manner with reference to formulae (3) and (4), but also 
by making a comparison between FIG. 3a representing 
g(k) of formula (3) and FIG. 3* representing g(k) in 
accordance with formula (4) for ko=2. Using the win- 
dow functions g(k) as defined in the formulae (3) and (4) 
consequently results in the same final value j3 of param- 
eter /3(ra) when the adaptation factor a has the same 
value for both cases. 

The convergence behaviour of the FDAF according 
to the invention can again be illustrated on the basis of 
simulation results for the same case as used to illustrate 
the convergence behaviour of the known FDAF's. 
Mamtaining the afore-mentioned stylizing, these simula- 
tion results are also represented in FIG. 4 by curve d, 
which relates to the FDAF including window means 11 
for implementing a window function g(k) in accordance 
with formula (4) with ko=0, the adaptation factor hav- 
ing a value a=2~ 5 for obtaining the prescribed final 
value /?= — 18 dB. This curve d substantially coincides 
with curve b in a large range from m=0 onwards and 
has the same final value j3 as curve b. Since curve b 
relates to the known FDAF with window means 11 
implementing a window function g(k) as defined in 
formula (3) and also having an adaptation factor 
a=2~ 5 , FIG. 4 is a clear illustration of the fact that the 
convergence behaviour of the FDAF according to the 
invention is in practice fully comparable to the conver- 
gence behaviour of this known FDAF. The deviations, 
shown in a stylized manner in FIG. 4, between curves b 
and d for low values /3(m) are illustrative of the fact that 
in the final stage the convergence proceeds somewhat 
slower for curve d than for curve b. This is however of 
little practical importance, as the convergence in the 
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final stage does indeed proceed slowly and the final 
value p is the same for both curves b and d. 

The preceding description can now be summarized as 
follows. For a FDAF in accordance with the overlap- 
save method, the article by Clark et aL describes a pre- 5 
ferred implementation containing five 2N-point FFTs, 
two of which are used to constrain the last N-time- 
domain weighting factors to be zero by using an appro- 
priate window function in the time-domain. By utilizing JQ 
a priori information about the impulse response to be 
modelled, the invention provides a FDAF in accor- 
dance with the overlap-save method having a special 
time-domain window function which can be imple- 
mented very efficiently in the frequency-domain, the 15 
FDAF only containing three 2N-point FFT's, whereas 
its convergence properties are comparable to those of 
the prior art implementation containing five 2N-point 
FFT's. 

What b claimed is: 20 
L A frequency-domain block-adaptive digial filter 
having a finite impulse response of length N for filtering 
a time-domain input signal in accordance with the over- 
lap-save method, comprising: 
first sectioning means for segmenting the time- 25 
domain input signal into blocks of length 2N, each 
block overlapping its preceding block over a 
length N; 

first transformation means for obtaining 2N frequen- 
cy-domain components of the Discrete Fourier 30 
Transform of length 2N of each input signal block; 

first multiplier means for multiplying each of the 2N 
frequency-domain components with an associated 
frequency-domain weighting factor for obtaining 35 
2N weighted frequency-domain components; 

second transformation means for obtaining 2N time- 
domain components of the inverse Discrete Fou- 
rier Transform of length 2N of the weighted fre- 
quency-domain components; 40 

second sectioning means for discarding the first N 
time-domain components and conveying the last N 



time-domain components as the time-domain out- 
put signal of the filter; 

means for generating a time-domain error signal as a 
difference between the filter output signal and a 
reference signal; 

third sectioning means for segmenting the time- 
domain error signal into blocks of length 2N, each 
block overlapping its preceding block over a 
length N, and for assigning the zero value to the 
first portion of a length N of each block; 

third transformation means for obtaining 2N frequen- 
cy-domain components of the Discrete Fourier 
Transform of length 2N of each error signal block; 

conjugation means for forming the complex conju- 
gate value of the 2N frequency-domain compo- 
nents produced by the first transformation means; . 

second multiplier means for multiplying each of the 
2N frequency-domain components produced by 
the third transformation means with the associated 
frequency-domain component produced by the 
conjugation means and with a gain factor for ob- 
taining 2N frequency-domain products; 

window means for performing an operation on the 
2N frequency-domain products, whose time- 

. domain equivalent is a multiplication by a window 
function, to obtain 2N frequency-domain 
weighting factor modifications; 

accumulator means for accumulating on a block-by- 
block basis each of the 2N weighting factor modifi- 
cations for obtaining the 2N frequency-domain 
weighting factors; 
characterized in that the window means are arranged 
for convolving the 2N frequency-domain products by a 
function having one real and two mutually conjugate 
complex coefficients and corresponding to a time- 
domain window function g(k) of length 2N defined by 

^fc)=a)[l + COS((Jk-*a)v/N] 

1 

for k=0, 1, . . . , 2N— 1, where ko is a constant with 
0Sk*<N. 
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